Resonant activation in a nonadiabatically driven optical lattice 
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We demonstrate the phenomenon of resonant activation in a non-adiabatically driven dissipative 
optical lattice with broken time-symmetry. The resonant activation results in a resonance as function 
of the driving frequency in the current of atoms through the periodic potential. We demonstrate 
that the resonance is produced by the interplay between deterministic driving and fluctuations, and 
we also show that by changing the frequency of the driving it is possible to control the direction of 
the diffusion. 

PACS numbers: 05.45.-a, 42.65.Es, 32.80.Pj 



o 



03 



i 

C 

O 
o 



> 

(N 

in 
in 

(N 

in 
o 

c3 



-a 
c 

o 
o 



X 



The problem of the escape of a Brownian particle out 
of a potential well, first characterized by Kramers pj, 
plays a central role in the description of many processes 
in physics, chemistry and biology [2j,|3j. Kramers' law 
predicts an escape rate of the form exp (—U/kT), where 
U is the depth of the well, T the temperature and k 
the Boltzmann constant. In the case of a nonadiabat- 
ically driven Brownian particle the aforementioned sce- 
nario may change significantly, and much work has been 
devoted to the study of this non-equilibrium phenomenon 
H H H 0, H H Efl EH . It has been shown that in the 
presence of nonadiabatic driving the lifetime of the par- 
ticle in the potential well can be significantly reduced, a 
phenomenon named resonant activation. 

Resonant activation was firstly observed in a current- 
biased Josephson tunnel junction Q, and more recently 
for a Brownian particle optically trapped in a double well 
potential [ljj. Resonant activation has also been theo- 
retically studied for Brownian particles in periodic po- 
tentials, a configuration relevant for the realization of 
Brownian motors [12j . Also in this case the nonadiabatic 
driving may result in a significant enhancement of the 
activation rate. Furthermore it has been predicted that, 
whenever the spatio-temporal symmetry of the system is 
broken, the resonant activation gives rise to resonant rec- 
tification of fluctuations 0, EH . Indeed, very recently 
Brownian motors have been realized with nonadiabati- 
cally driven Brownian particles [T3 . [l5| . However in those 
works the resonant nature of the rectification mechanism, 
and the underlying resonant activation process, have not 
been investigated. In particular, the possibility to control 
the amplitude and direction of the diffusion by changing 
the frequency of the driving field, as predicted by theo- 
retical work j(| 0, ^3 > nas n °t yet been demonstrated. 

In this Letter we demonstrate experimentally the phe- 
nomenon of resonant activation in a Brownian motor by 
using cold atoms in a driven dissipative optical lattice as 
a model system. Indeed, we observe the appearance of a 
resonance while monitoring the current of atoms through 
the lattice as a function of the driving frequency. We 
show that this resonance is due to the interplay between 
fluctuations and deterministic driving, and we demon- 
strate that by varying the driving frequency it is possible 



to reverse the current direction. 

In our experiment we use cold atoms in a dissipative 
optical lattice |l6l |. in which the atom- light interaction 
determines both the periodic potential for the atoms and 
the dissipation mechanism which leads to a friction force 
- the so called Sisyphus cooling - and to fluctuations in 
the atomic dynamics. This system offers the significant 
advantage of being easily tunable over a wide range of pa- 
rameters: the potential depth, the fluctuations level and 
the parameters (frequency and amplitude) of the driving 
force can be varied and carefully controlled over a large 
interval of values. This is an essential feature for the 
investigation presented here. 

Before presenting the experimental results, we ana- 
lyze our system through numerical simulations. For 
the sake of clarity, our theoretical analysis will be lim- 
ited to the simplest one-dimensional configuration in 
which Sisyphus cooling takes place: a J g — 1/2 — > 
J e = 3/2 atomic transition and two counterpropagat- 
ing laser fields with orthogonal linear polarizations - 
the so called lin_Llin configuration. The light interfer- 
ence pattern results in two optical potentials U± for the 
atoms, one for each ground state |±) , in phase opposition: 
U± = Uq [— 2 ± cos 2kz] /2, where z is the atomic posi- 
tion along the axis Oz of light propagation, k the laser 
field wavevector and Uo the depth of the optical poten- 
tial. The stochastic process of optical pumping transfers 
the atom from one ground state to the other one, chang- 
ing in this way the optical potential experienced by the 
atom. This stochastic process results in a friction force 
and produces fluctuations in the atomic dynamics |16| . 
The departure rates 7±^=f(z) from the |±) states can 
be written in terms of the photon scattering rate T' as 
7±^ T = r'(l ±cos2fcz)/9 |l6(. It appears therefore that 
the amplitude of the fluctuations can be quantitatively 
characterized by the photon scattering rate T', which is 
an experimentally accessible parameter. 

To study the phenomenon of resonant activation we 
need to drive the atoms with a zero-mean oscillating 
force. We consider here an ac force consisting of two 
harmonics, A\ cos (ujt) and cos {2ojt — tt/2), so that 
the resonant activation should lead, following the break- 
ing of the time-symmetry of the system, to a resonant 
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generation of a current, as predicted by theoretical work 
HI H> • I n the numerical work it is obviously possi- 
ble to "apply" directly an homogeneous ac force to the 
atoms, by simply including the appropriate terms in the 
equation of motion. On the contrary, in the experiment 
this is not possible, and forces can be applied only by 
phase modulating the lattice beams. For consistency, we 
follow the same approach in the theoretical analysis, and 
we consider a phase-modulation a(t) of one of the lattice 
beams of the form 



a(t) = a Q 



Ai cos (tut) H — A2 cos (2uit — 7r/2) 



(1) 



In this way in the accelerated frame in which the optical 
potential is stationary the atoms experience an inertial 
force 

2 

F(t) = muj a ° M 1 C o S f w t) + A 2 cos (2wt - tt/2)] , (2) 

where m is the atomic mass and k the laser field wavevec- 
tor. 

To study the atomic dynamics in the presence of the 
nonadiabatic driving, we follow the same procedure de- 
veloped to investigate laser cooling processes in (un- 
driven) optical lattices. The Fokker-Plank-type equation 
for the undriven system, and the Monte Carlo technique 
to derive the atomic trajectories have been discussed 
in detail in Refs. (liil Hq . The generalization of that 
method for the driven system of interest here is straigh- 
forward, and consists in the inclusion of a time-dependent 
shift a(t) (see Eq. (QJ) in the relative phase between the 
two laser fields generating the optical lattice. 

Through Monte Carlo simulations, we calculated the 
mean atomic velocity v as a function of the frequency ui 
of the driving, for different amplitudes of the ac force. 
The results of our calculations are shown in Fig. ^ For 
each data set, in order to keep constant the amplitude 
Fq = muj 2 ao/2k of the ac force (see Eq. we varied 
the amplitude «o of the phase modulation according to 
aa = a/uj 2 , with a constant for a given Fq. This is the 
same procedure used in the experiment. 

Figure ^ shows clearly the appearance of a resonance 
in the current amplitude as a function of the driving fre- 
quency. Indeed, for weak driving the current is negligible. 
At larger amplitude of the driving the current differs sig- 
nificantly from zero, and shows a well defined resonance. 
The resonance is observed in the regime of non-adiabatic 
driving, i.e. for driving frequencies of the order of or 
exceeding the vibrational frequency. The numerical sim- 
ulations also show that by changing the driving frequency 
it is possible to reverse the current direction, as predicted 
by the general theory 0, 0, [n| . We note that we carried 
out numerical simulations for two different ratios of the 
force harmonics' amplitude: A1/A2 = 1/4 (Fig. and 
A1/A2 = 1 (not shown). The two sets of calculations 
evidenced the same qualitative behaviour. 

To determine the nature of the resonance we focus our 
attention on the range of frequencies of the ac fields cor- 
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FIG. 1: Results of Monte Carlo simulations for a sample of 
n = 10 4 atoms in a ID lin_Llin optical lattice. The mean 
atomic velocity v, rescaled by the atomic recoil velocity v r = 
hk/m, is shown as a function of the driving frequency ui, for 
different amplitudes of the ac force. The driving frequency u 
is expressed in terms of the vibrational frequency U) v of the 
atoms at the bottom of the potential wells. Parameters of the 
calculation are: T' — 10u> r , Uo ~ 100 ■ E r , with E r and ui r 
the recoil energy and frequency, respectively. The amplitude 
of the two harmonics of the force are: A% = 1, A% = 4. 



responding to non-adiabatic driving. This regime is il- 
lustrated in Fig. [2 where the resonance in the atomic 
current as a function of the driving frequency is shown 
for a given amplitude of the ac field and at different val- 
ues of the scattering rate. From Fig. |2 it appears that 
the amplitude of the resonance shows a non-monotonic 
dependence on the scattering rate T': at small V the 
resonance amplitude increases with V , then reaches a 
maximum and finally at larger values of T' decreases. 
The constructive role played by the noise at low levels 
of T' shows that the resonance observed in the numer- 
ical simulations is determined by the interplay between 
the applied ac forces and the random fluctuations, which 
results in the rectification of the latter ones. This is at 
variance with the behavior observed at frequencies some- 
what below uj v where the magnitude of the (negative) 
current is decreased by an increase of the scattering rate, 
a behaviour characteristic of deterministic rectification of 
forces. The conclusion of our numerical analysis is there- 
fore that the resonant activation phenomenon should be 
observable in a dissipative optical lattice, and should re- 
sult in a resonance in the atomic current as a function of 
the driving frequency. 

The experiment is a direct implementation of what de- 
scribed in the theoretical analysis. Instead of using a ID 
optical lattice, as in the numerical work, we use a 3D 
lattice, which offers the significant advantage to confine 
the atoms in the three directions. This reduces the losses 
of atoms from the lattice during the experiment. 
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FIG. 2: Results of Monte Carlo simulations for a sample of 
n = 10 4 atoms in a ID lin_Llin optical lattice. The mean 
atomic velocity v, is shown as a function of the driving fre- 
quency iv, for different values of the scattering rate V . In 
the inset, the mean atomic velocity u P cak at the maximum of 
the resonance is reported as a function of the scattering rate. 
The depth of the optical potential is Uo = 100 ■ E r ; the ac 
force amplitude corresponds to a — IQQ-klu^.. The amplitudes 
of the two harmonics of the force are : A\ = 1, Ai — 4. The 
lines are guides for the eye. 



The experimental set-up is the same as the one used 
in Ref. |15|. and consists of four linearly-polarized laser 
beams arranged in the so-called umbrellalike configura- 
tion [l^. One laser beam propagates in the z-direction. 
This is the beam which will be phase-modulated. The 
three other laser beams propagate in the opposite di- 
rection, and are arranged along the edges of a triangu- 
lar pyramid having the z-direction as axis, with the az- 
imuthal angle between each pair of beams equal to 27r/3. 
For this lattice beam configuration, the interference of 
the laser fields produces a periodic and spatially sym- 
metric optical potential, with the potential minima asso- 
ciated with pure circular (a + or er~) polarization of the 
light [l8|]. For an atom with a, F g = F ^ F e = F+l tran- 
sition, the optical potential consists precisely of 2F + 1 
potentials, one for each ground state sublevel of the atom. 

Cesium atoms are cooled and trapped in a magneto- 
optical trap (MOT). At a given instant the MOT is 
turned off and the four lattice beams are turned on. 
The atoms are left in this undriven optical lattice for 
2 ms. This is sufficient for the atoms to thermalize and 
reach equilibrium. Then the phase modulation a(t) [see 



Eq. Q] is slowly turned on. The dynamics of the atoms 
is studied with a charged-coupled device (CCD) camera. 
After a short transient, the center-of-mass of the atomic 
cloud is observed to be set into uniform motion along the 
z axis, and a center-of-mass velocity v is correspondingly 
derived. 
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FIG. 3: Experimental results for the average atomic velocity 
as a function of the driving frequency, for different ampli- 
tudes of the driving force. The optical potential is the same 
for all measurements and corresponds to a vibrational fre- 
quency U) v — 2n ■ 170 kHz. The driving frequency satisfying 
the condition 2ui = uj v is indicated by an arrow. The detuning 
A of the lattice from atomic resonance is A = 11. IT, where 
r is the excited state linewidth. The values for the velocity 
are expressed in terms of the recoil velocity v r , equal to 3.52 
mm/s for the Cs D2 line. The two harmonics of the force have 
equal amplitude: A\ = A2 — 1. The lines are guides for the 
eye. 

Results for the average atomic velocity as a function 
of the driving frequency are shown in Fig. [21 for different 
values of the amplitude of the driving. We clearly observe 
the build-up of a resonance when the amplitude of the 
driving is progressively increased. The resonance appears 
in the regime of non-adiabatic driving (2a; >u„), and a 
current reversal is observed on the low-frequency side of 
the resonance, in agreement with our simulations and 
with the general theory ||| || |n| . 

We examine now the dependence of the resonance on 
the scattering rate, i.e. the fluctuations amplitude, with 
the aim to demonstrate that the observed resonance is 
indeed produced by the interplay of deterministic driv- 
ing and fluctuations. Figure 21 shows our results for 
the average atomic velocity as a function of the driv- 
ing frequency, at different values of the scattering rate. 
It clearly appears that the resonance amplitude shows a 
non-monotonic dependence on the scattering rate: the 
resonance initially increases at increasing values of T', 
then reaches a maximum and starts decreasing at large 
values of the scattering rate. This is in agreement with 
our numerical results and clearly demonstrates that the 
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FIG. 4: Experimental results for the average atomic veloc- 
ity as a function of the driving frequency, at different scat- 
tering rates. In the inset the resonance amplitude, i.e. the 
peak mean atomic velocity, is reported as a function of the 
scattering rate. The optical potential, and the driving force 
amplitude, are the same for all data sets, and they are char- 
acterized by w v = 2tt • 170 kHz and a = 13.3(MHz) 2 . The 
driving frequency satisfying the condition 2u = uj v is indi- 
cated by an arrow. Different data sets correspond to different 
detunings A, i.e., to different scattering rates as the optical 
potential is kept constant. The data are labeled by the quan- 
tity r s = [cji,/(27r)] 2 /A which is proportional to the optical 
scattering rate. The two harmonics of the force have equal 
amplitudes: A\ = A 2 = 1. The lines are guides for the eye. 



resonance is determined by the interplay between deter- 
ministic driving and fluctuations and, due to the broken 
time-symmetry, results in the rectification of the latter 
ones. 

In conclusion, in this work we demonstrated the 
phenomenon of resonant activation in a nonadiabati- 
cally driven dissipative optical lattice with broken time- 
symmetry. Due to the broken symmetry of the system, 
the resonant activation results in a resonance in the cur- 
rent of atoms through the periodic potential. We demon- 
strated that the resonance is produced by the interplay 
between deterministic driving and fluctuations, and we 
also showed that by changing the frequency of the driv- 
ing it is possible to control the direction of the diffusion, 
as predicted by theoretical models. We notice that the 
rectification of fluctuations with non-adiabatically driven 
Brownian particles was already observed in previous work 
[l4lll5j . but the resonant nature of the rectification mech- 
anism was not demonstrated. Our work therefore also 
establishes experimentally the connection between reso- 
nant activation and resonant rectification of fluctuations, 
and confirms the theoretical predictions. The present ex- 
perimental realization, in which both deterministic and 
fluctuating forces originate from light fields, also shows 
the generality of the phenomenon of resonant activation, 
which is not restricted to systems in which the fluctua- 
tions are of thermal origin, as usually considered in the- 
oretical work. 

We thank EPSRC, UK and the Royal Society for fi- 
nancial support. 



[1] H.A. Kramers, Physica 7, 284 (1940). 

[2] P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 
62, 251 (1990). 

[3] V.I. Melnikov, Phys. Rep. 209, 2 (1991). 

[4] M.H. Devoret, J.M. Martinis, D. Esteve, and J. Clarke, 
Phys. Rev. Lett. 53, 1260 (1984); M.H. Devoret, D. Es- 
teve, J.M. Martinis, A. Cleland, and J. Clarke, Phys. 
Rev. B 36, 58 (1987). 

[5] R. Mannella and P. Grigolini, Phys. Rev. B 39, 4722 
(1989). 

[6] M.I. Dykman, H. Rabitz, V.N. Smelyanskiy, and 
B.E. Vugmeister, Phys. Rev. Lett. 79, 1178 (1997). 

[7] D.G. Luchinsky et al, J. Phys. A: Math. Gen. 32, L321 
(1999). 

[8] D.G. Luchinsky, M.J. Greenall, P.V.E. McClintock, 

Phys. Lett. A 273, 316 (2000). 
[9] J. Lehmann, P. Reimann and P. Hanggi, Phys. Rev. Lett. 

84, 1639 (2000); Phys. Rev. E 62, 6282 (2000). 
[10] M.I. Dykman et al, Chaos 1, 587 (2001). 
[11] S.M. Soskin, R. Mannella, and P.V.E. McClintock, Phys. 

Rep. 373, 247 (2003). 
[12] For early investigations on Brownian motors, see: 

M.O. Magnasco, Phys. Rev. Lett. 71, 1477 (1993). A. Ad- 

jari et al, J. Phys. I (France) 4, 1551 (1994); M.C. Ma- 



hato and A.M. Jayannavar, Phys. Lett. A 209, 21 (1995); 
D.R. Chialvo and M.M. Millonnas, ibid. 209, 26 (1995); 
R. Bartussek, P. Hanggi and J.G. Kissner, Europhys. 
Lett. 28, 459 (1994). For a review on Brownian mo- 
tors see R.D. Astumian and P. Hanggi, Phys. Today 
55, 33 (2002); P. Reimann, Phys. Rep. 361, 57 (2002); 
P. Reimann and P. Hanggi, Appl. Phys. A 75, 169 (2002); 
see also P.V.E. McClintock, Nature 401, 23 (1999). 
[13] I. Goychuk and P. Hanggi, Europhys. Lett. 43, 503 
(1998). 

[14] S. Weiss, D. Koelle, J. Miiller, R. Gross, and K. Barthel, 

Europhys. Lett. 51, 499 (2000). 
[15] P.H. Jones, M. Goonasekera, and F. Renzoni, Phys. Rev. 

Lett. 93, 073904 (2004). 
[16] For a recent review of optical lattices, see G. Grynberg 

and C. Mennerat-Robilliard, Phys. Rep. 355, 335 (2001). 
[17] For a general analysis of the role of the symmetries in 

the case of bi-harmonic driving, see: S. Flach, O. Yev- 

tushenko and Y. Zolotaryuk, Phys. Rev. Lett. 84, 2358 

(2000); S. Denisov et al, Phys. Rev. E 66, 041104 (2002). 
[18] K.I. Petsas, A.B. Coates, and G. Grynberg, Phys. Rev. 

A 50, 5173 (1994). 



